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Abstract 

Starting from a relativistic s-wave scattering length model for the two par- 
ticle input we construct an unambiguous, unitary solution of the relativistic 
three body problem given only the masses m a ,mb,m c and the masses of the 
two body bound states fJ-bc, f^ca, Hab- 
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1 INTRODUCTION 



Some time ago one of us made an analysis of what was then needed to "solve" the 



three nucleon problem ||17||. The conclusions were rather discouraging. In particular, it 
appeared to us that the enormous effort that Tjon and his collaborators, and a number 
of other few body nuclear physicists, were making to obtain and use "realistic" two- 
nucleon potential models to calculate the triton wave function and related problems 
were necessarily plagued by ambiguities which it would be difficult to remove. We 
eventually came to the conclusion that a "zero-range" or "on-shell" treatment of the 
problem might aid in simplifying the contact between these efforts and experiment. 
Unfortunately our efforts did not lead to the desired result in the form of a unique 
theoretical analysis for reasons we discuss to some extent in our paper reporting 



this body of research^]. The reasons are rather complicated, but can be briefly 
summarized by the statement that only a finite particle number relativistic scattering 
theory can provide the needed framework. This paper presents a relativistic scattering 
length model which we believe will be a practical first step toward creating such a 
fundamental theory. 

The key to understanding how we can create a simple, soluble and unitary model 
for a relativistic three body system is to realize that we embed the two body subsys- 
tems in the three body space from the start. We consider two particles with invariant 
masses m a , m b which scatter elastically when Mandelstam variable s (square of the 
invariant four momentum) lies in the range (m a + m b ) 2 < s < M 2 h ] here M t ^ is the 
total energy of the first inelastic (particle creation) threshold in the zero 3-momentum 
reference frame and k(s) is the magnitude of the momentum of either particle in this 
frame. On 3-momentum and energy shell, we write s = M 2 with the conventional al- 
gebraic connections to the energies of the free particles outside the elastic scattering 
volume: 

2 2 2 2 2 

e a~ m a = P =e b -m b 

M 2 = (e a + e b ) 2 - \p a +p b \ 2 

m M„) = [(^-("■.whap -(*■.- *.)')]» (1) 
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(m a + m b ) 2 < s = M 2 < M? t 



th 



As in our non-relativistic treatment |21[ , we consider only an interaction in which 
no angular momentum is transferred. We also analytically continue to values of 
s 7^ M 2 , keeping 3-momentum conserved as in the non-relativistic case, but using this 
kinematics and two-particle amplitudes embedded in a multi-particle space. As was 
pointed out to one of us (HPN) by J.V.Lindesay 0] in commenting on a preliminary 



version of this paper, the correct way to carry out this analytic continuation is to 
take 



. L „ - (m a + m b ) 2 ] [s - (m a - m b ) 2 } 
k ab {s;m a + m b ) = +]l — (2) 

This has the effect of embedding our two-particle interaction in a multi-particle space. 
We were earlier led to a version of this "off shell extension" by other considerations fl6|| . 
We are much indebted to Castillejo0 for pointing out to us in a discussion of [[LJJ 
that we had, in fact, implicitly assumed we were in a multi-particle space. 

Whatever the off-shell extension, we can insure on-shell unitarity for the scattering 
amplitude T(s) with the normalization Im T(s) = £;(s)|T| 2 in the elastic scattering 
region by using the scattering length formula 

_ 1 _ e^sin 5(s) 
I[S) ~ 7 - ik(s) ~ k(s) [6) 

where 7 is any finite constant. To generate a bound state pole at s — [i 2 b we can (for 
H 2 ab > (m a - m b ) 2 ) take 

4/47 2 = [(m b + m c f - //J \^ ah - (m a - m b ) 2 } > (4) 

and pick the branch in the square root defining k(s) by analytic continuation below 
elastic scattering threshold to insure that a pole occurs p5| . As we discuss to some 



extent in pl| , p. 1869 and as Weinberg explored extensively in his quasi-particle ap- 



proach to the three body problem |29|, |30| , plj , this prescription amounts to assuming 
that the two particles composing the bound state are "structureless" . This completes 
our specification of the two particle input to our three body problem. 

In the next section we show that in our context this algebraic construction of on- 
shell two body unitarity plus the algebraic form of the Faddeev equations guarantees 
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the three body unitarity of their solution. This was not obvious to some people in 
the context of our previous non-relativistic analysis, which led to considerable formal 



complexity in the final presentation In the current context the triviality becomes 
manifest for the 2-2 channels below three body breakup threshold where the solution 
is algebraic, as we show in Sec. 3.4. The 3 free to 3 free, coalesence and breakup 
amplitudes require a more detailed analysis of the zero angular momentum equations 
than we initially intended to make here. However, in order to meet a referee comment 
on an earlier version of this paper f23|, we have extended our discussion in Sec. 3 and 
there provide algebraic solutions of the separable integral equations which follow 
from our assumptions when the three free particle channels are explicitly included. 



Examples of these solutions will be explored elsewhere |fL5|| . Our concluding section 
contains speculations as to where this simple result might be applied. 



2 RELATIVISTIC THREE BODY UNITARITY 

2.1 Relativistic Faddeev 2-2 channel kinematics 

Our model contains three structureless particles labeled a, b, c, which taken pairwise 
can form three bound states (be) , (ca) , (ab) whose only structure arises from these 
simple constituents. Although we are dealing with a relativistic system, we restrict 
our energies to the range which allows no particle creation; in this respect the situation 
is the same as in the non-relativistic three body problem, and the Faddeev[|. |J 
channel decomposition can be employed. Because of this simple structure, if we start 
in a state with zero total 3-momentum and angular momentum, this situation will 
persist so long as our interactions have no internal degrees of freedom. This is true by 
hypothesis for the model we described in the introduction. We further assume that 
we start the system out with a scattering between one of the particles (which we can 
pick to be a) and the implied bound state pair (which will then be be). Following 
Faddeev, we drop the redundant index. We distinguish the bare particles from the 
bound states by using the symbols m and \x for their respective masses. Then a 
specific scattering problem is completely specified by supplying numerical values for 
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the following eight parameters: 

Particle masses : m a , m b , m c 

Bound state masses : u a , u b , u c (5) 

Invariant 4 — momentum : M = 

Input momentum : k a 

The input momentum is the magnitude of the 3-momentum of particle a of mass m a in 
the zero momentum frame with the bound state of particles b and c of mass u a having 
3-momentum of equal magnitude but opposite direction. We call this momentum n a 
Calling the corresponding energy e a and the energy of the bare particle e a fixes the 
exterior kinematics as follows: 

e l ~ k l = m li xea,b,c 

4 - *4 = & x ea,b,c (6) 

M abc = e x + e x = M independent of x 

As can be seen from Eq.7 below, the constraint M > m x + u x is required, introducing 
the ^-functions explicitly noted in Eq.52 below which mark the opening of the 2+1 
channels at their energetic thresholds; these ^-functions are required for consistency 
with the orthogonality (Eq.45) and completeness (Eq.46) relations needed in the most 
general case we discuss. 

We will not need to leave the three body zero momentum frame in this paper. 
Given M, we can in this restricted environment immediately compute the input mo- 
menta because 

4M 2 k 2 x (^ x ,M)= [M 2 -(m x + u x ) 2 ][M 2 -(m x -fi x ) 2 ] =4M 2 K 2 X (7) 

x G a,b,c 

Thus the input momentum is not a parameter, but we still have to specify whether the 
input channel is a or b or c. Above three body breakup threshold (M > m a + m b + m c ) 
we must specify both the input and the output momenta for the "spectator" (see Sec. 
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3.1). Because we can have one or two rearrangement output channels in addition to 
elastic scattering in the entrance channel, the entrance channel has to be open for 
anything to happen, but how many output channels will then be reached introduces 
an asymmetry into the predictions. The next problem is to insure that these various 
possibilities are described by a unitary formalism. 



2.2 On Shell Faddeev Equations 

If we start in the a channel, the simplest thing that can happen is that the be pair 
scatter without the initial spectator a being affected. These are the "disconnected 
diagrams" which initially caused so much trouble in analyzing the three body scat- 
tering problem, until Faddeev realized that that they could be subtracted out of the 
amplitude, leaving only connected diagrams. Here we call the full amplitude 



(3) 



T(M) = Z x>yeaM ®T xy {M)9{M -m x - fi x )9(M - m y - fx y ) 



where {3) T xy (M) are the unique solutions of Eq. 10 below. 

We will return to the functions in Sec. 3, where we compute only those pro- 
cesses allowed by the conservation laws, except when a bound state vertex opens or 
closes in a single channel. We call the disconnected amplitude ^T a (M) and define it 
by 

®T a (M) = ^T bc {[M - e a {M)f) (9) 

where e a (M) is to be computed from Eq. 6. Then the relativistic version of the 
"on-shell Faddeev Equations" for our simple pole model become 



Tab — T a 5 a b = T a RTbb + T a RT c b = T a T, x 5a X RT x b 

= T aa RT h + T ac RT b = \L y T ay R5 yh ]T h (10) 
cyclic on a, b, c 

£>ab = 1 — 8 a b 

Here R is the three free particle propagator, which because of our on-shell kinematics 
is simply a constant whose value we will determine from unitarity in the next two 
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sections. The fact that the two alternative forms of Eq. 10 define the same function is 



critical for time reversal invariance. That requirement was what led us to conclude[21 



that our non-relativistic zero range model could not be consistently defined when the 
two body amplitudes contain "left hand cuts" and not just bound state or CDD[^] 
poles. Algebraically the fact that summing the multiple scattering series starting 
from either the first or the last scattering — which is what the two alternative forms 
of Eq.10 express — follows from the fact that the two forms define the same multiple 
scattering series. Convergence follows if we use two-body input amplitudes which are 
always less than unity in absolute magnitude [ITjJ]. 

In our theory the convergence of the two forms of the multiple scattering series 
given in Eq. 10 to the same function suffices to prove the uniqueness of the solution. 
In a normal Hamiltonian theory, one would require that the homogeneous solutions 
of the two equations also must be proved to be identical. A related question is to ask 
"what operator must be diagonalized in order to make the spectral expansion (46)." 
Both questions are discussed more fully in the non-relativistic context when taking 
the "zero range limit" (see Ref. 21). Note that the same dispersion relation used to 
establish the unitarity of the two-body input in the non-relativistic case also applies 
to the relativistic scattering length model used here. 

There is no known unambiguous way to go from any field-theoretic description of 
strongly interacting particles to a non-relativistic potential model of the class used 



phenomenologically in nuclear physics [|17]]. Ref.21 tried to remedy this defect by start- 
ing from a dispersion- theoretic representation of the two body amplitudes which is 
the non-relativistic limit of dispersion-theoretic relativistic S-Matrix amplitudes (cf. 
Ref.25), leading to a "left hand cut" representing virtual relativistic particle exchanges 
(eg pions). The difficulty which showed up in implementing this idea was that this 
approach cannot be made general enough. In particular, the Low equation (which, 
superficially, would seem to allow us to construct the non-relativistic potential from 
knowledge of the two-body phase shifts and the left hand cuts in the correspond- 
ing amplitudes) cannot be constructed. Specifically, the requirement of time reversal 
invariance for the two-body potential model which results from the attempted con- 
struction cannot be satisfied when the dispersion-theoretic input has left hand cuts. 
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Consequently any bound state pole in the two body-input will not vanish in the zero 
range limit, and corresponds to a CDD pole (Ref.7) as noted above. Consequently, 
some relativistic theory must be invoked to provide the parameters of this pole. As 
mentioned in Sec. 3. 4 we use here the "handy-dandy formula" connecting masses to 
coupling constants. We conclude that there is no direct way to relate our model to a 
non-relativistic Hamiltonian model, and hence that our completeness relation is part 
our model; it need not be derived. 

The second point to emphasize is that our two-body bound states are structureless. 
They only appear in our three body system above energetic threshold, when the 
third particle can play the kinematic role of a "spectator". For the same reason, 
the equations in this paper cannot lead directly to three body bound states, which 
is why we do not need the (non-existent) "homogeneous solutions" to establish the 
equivalence of the two orders of the operators in Eq.10. The convergence of the 
summation of the two multiple scattering series to the same result suffices. Consistent 
treatment of three body bound states requires us to embed this model in a four body 
space in order to provide a kinematic "spectator" , and is not attempted in this paper. 



2.3 Three Particle Unitarity from Two Particle Unitarity 

That the unitarity of the off-shell Faddeev equations follows from the unitarity of the 
two body input amplitudes was shown by Freedman, Lovelace and Namyslowski flT 



and independently by Kowalski[13|, who taught HPN the simple algebraic proof given 
below. The key is to write the unitarity condition on the two-body amplitude t a in 
the three body space as 

t a (M) - t* a (M) = t a (M)(R - R*)t* a (M) (11) 

which determines the normalization of R in terms of the normalization of t a . In order 
to avoid kinematic factors, it is convenient to use 

t a (M) = [e i5 ^ s) sin 5 a (s)] s=[M -e a (M)]i = [hc(s) (2) T 6c (s)] s=[M _ ea w (12) 

so that t — t* = 2i tt* . Then R = +i provides a channel independent propagator 
in the three particle space. The proper thresholds for the opening of the various 
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channels are provided by the functions in Eq. 8. Using the same normalization, 
the unitarity condition on the three body channel amplitudes is then simply 



T(M) - T*(M) = T(M)(R - R*)T*(M) (13) 

and below breakup threshold the on-shell Faddeev equations become the algebraic 
equations 

Tax ta^ax ~\~H a (T bx ~\~ T cx ^j -\-it a Eiy5 a yTy X 

= +i((T ay + T az )t x = +i[T>yT ay 5 yx }tx (14) 



The proof of unitarity follows the same steps taken inpO|, Eq.'s (2)-(5), but below 



breakup threshold these are now actually algebraic rather than symbolic. If in T ab we 
call L a = SfcT a b and Ft, = T, a T a b, we clearly have that 

T(M) = S a L a = Y, b F b (15) 

Then the unitarity condition we wish to prove becomes 

rp p* ? 

Z X F X (R-R*)L* X + E x J xy F x (R-R*)L* y (16) 

where we have separated out the term where the indices differ so that we can make 
use of the two body unitarity condition. This is possible because, in the current 
notation, we can rewrite the Faddeev equations (Eq.10) as 

F x (1 ^->wF w S wx R)tx 

L* x = t* x (l-E z R% x L* z ) (17) 
Consequently the equation to be proved becomes 

rp p* ? 

£z(l — ^wF w 5 wx R) t x (R — R*)t* x (1 — T lz R*5 zx L* z ) (18) 

+ E>x,yfixyF x {R* — R)L* y 

9 



We can now take the critical step of substituting t x — t* for t x (R — R*)t* x and find 
that 

Sjb(1 — T lw F w 5 wx R)t x — t*(l — E Z R*S ZX L* Z ) (19) 

^wRw^wxRtx Fx\R $xyLy 
+ ^x,yF x RS xy [t* — t*E 2 i?*^a;L* — L*] 

T — T* Q.E.D. 

where the unwanted terms vanish because F x and L* are solutions of the Faddeev 
equations. 

3 SOLUTION OF THE ON-SHELL FADDEEV 
EQUATIONS 

3.1 3-3 Kinematics 

In order to properly relate the amplitudes for elastic and rearrangement collisions 
below and above breakup threshold it is useful to first express the amplitudes as 
operators in an orthonormal and complete space and then reduce them to integral 
or algebraic equations. In our non-relativistic treatment [pi]] this happened automati- 



cally because there we took the "zero range limit" of equations originally written in a 
larger space. Here we make a new approach by formulating the operators directly in 
the on-shell space of empirically observable particle and bound state momenta. For 
3-3 collisions, under the restriction to s-channel driving terms and total angular mo- 
mentum zero (and ignoring for the moment the two-body bound states and channels 
associated with them), this is simply the space of the relativistic Dalitz plot. We 
modify the standard notation to accommodate the Faddeev channel decomposition 
(a the spectator of a be pair interaction or scattering) as follows. For single free 
particles pi = (E^pt), % 6 a,b, c, and pf = Ef — Pi ■ Pi = mf. Then, for i, j, k cyclic 
or anti-cyclic on a, b, c, if Pi is the four-momentum of the spectator, we define the 
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Mandelstam invariant s$ for this channel by 

Si = m 2 k = (pj + p k ) 2 = (Ej + E k ) 2 - (fij + p k ) ■ (pj + p k ) (20) 

In an arbitrary coordinate system 

P = p a +p b + p c ; P 2 = M 2 = (E a + E b + E c ) 2 - (p a + p b + p c ) ■ (p a + p b + p c ) (21) 

It follows immediately that in the coordinate system where P = (p a + p b + p c ) = 0, 
that 

Si = M 2 + m 2 - 2MEi 
s a + s b + s c = M 2 + m 2 a + m 2 b + m\ = S 2 (22) 

and that, cf (34.20b), with \p a \ 2 -> fc 2 , 

AM 2 k 2 a (s a ) = [M 2 -(sl+m a ) 2 ][M 2 -(sl-m a ) 2 ] 

= (M 2 -m 2 a - s a f -Amls a (23) 
= (M 2 + m 2 a - s a ) 2 - 4M 2 m 2 a 

allowing us to define k a (s a ) = \p a \ in this zero 3-momentum frame. 

The magnitudes of the three momenta k a , k b , k c , or equivalently (cf. Eq.'s 22 and 
23) the three invariants s a ,s b ,s c , specify a rigid triangle. The orientation of this 
triangle relative to some space fixed system of axes in which the system as a whole 
has zero 3-momentum can be specified by three Euler angles, (a, (3, 7), cf. Eq. 34.18. 
In, for example, Osborn's treatment []27|] , the three degrees of freedom connecting 
the initial to the final state are discretized as a rotation matrix d J xx ,{9) and cos 9 
is expressed in terms of the scalar momenta k, k'. Under our s-channel, zero total 
angular momentum assumption, this rotation matrix is simply a constant, which is a 
way of seeing why we have only three degrees of freedom even before we go on shell 
(restrict ourselves to the interior of the Dalitz plot by fixing M in Eq.22). 
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3.2 Separable integral equations for 3 free- 3 free scattering 

Define the operators M. ax = %,8ax + Wax, making the Faddeev equations starting 
from entrance channel a read 

w aa - T a n(w ba + w ca ) = o 

-T b KW aa + W ba - T b KW ca = %KT a (24) 

-%n(w aa + w ba ) + w ca = r c TiT a 

Separate W ax into the four amplitude operators 

Wax A ax ~\~ $ax ~\~ C- ax ~\~ T^ax (^^) 

A ax '■ Anelastic scattering 

B ax : Breakup 

C ax '■ Coalesence 

V ax : 3 free — 3 free scattering 

Assume that the first three are zero. Then the operator Faddeev equations for the 
3-free to 3-free amplitude starting in entrance channel a are 

v aa - r a n(v ba + v ca ) = o 

-T b UV aa + V ba -T h nV ca = T b UT a (26) 
-%K(V aa + V ba ) + V ca = T c UT a 

and the two obvious cyclic permutations for entrance channels b and c. The alternative 
form, with exit channel (rather than entrance channel) specified as a is 

V aa -{V ab + V ac )TZT a = 
-V aa 1lT b + V ab -V ac KT b = T a TZT b (27) 
-(V aa + V ab )KT c + V ac = T a UT c 

To convert these into one-variable integral equations, we assume the orthogonality 
and completeness relations 

<s x \s' y > = 8 xy S(s x - s' y ) (28) 



I, 



(M-m x f 

x I ^X *^* > < ~^~ ^X I 

(m y +m z ) 2 
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for the states onto which we project. Here, and from now on, the kinematic constraint 
(cf. Eq.22) on the variables 

s a + s b + s c = M 2 + m 2 a + m\ + m\ = T? = s' a + s' b + s' c (29) 

is to be understood. The matrix elements of the operators are taken to be 

< s x \1Z\s' y > = R5 xy 9(M - m a - m b - m c ) 

< s x \T x \s' y > = 5 xy 5(s x - s' y )e l5x sin 5 x (s x , M)6(s x - (m y + m z ) 2 ) (30) 
< s x \V xy \s' y > = D xy (s x ; M; s' y )6(M - m a - m h - m c ) 

where R is a constant to be fixed by requiring on-shell unitarity and 5 xy — 1 — 
5 xy . It is convenient to also define t x (s x ,M) = e l5x sin S x (s x ,M). Then, by taking 
matrix elements of the operator equations between the appropriate bras and kets, 
and invoking the completeness relation to bring in the appropriate integrals, we find 
that, for a the entrance channel, 

D aa (s a , s' a ) - t a (s a )R[J^ a - + Z% ds'lD ba {sl s' a ) + jfcZ$ 2 ds»D ca (s'>; <)] 

= 

-t b (s b )RjlZ~ + Zi ds'paaK; <) + D ha (s b ; <) - t b (s b )Rf^- + Z$ 2 ds^D ca {s^ s' a ) 

= t b (s b )Rt a (s' a ) (31) 
-UsMCZi ds><D aa {sl s' a ) + jIZTS ds'lD ba {sl s' a )\ + V ca {s c] s' a ) 

= t c (s c )Rt a (s' a ) 

Here we have suppressed the fixed argument M in the functions and retained only 
the variables s x , s' x , s" of the integral equation and the fixed parameter s' a describing 
the initial state. The alternative form of the equations (Eq. 27) reverses the role of 
variables and parameters, leading to 

D aa (s> a ; s a ) - [SZrS ds'lD ab {s> a - s'l) + ds"D ac (s' a ; s^]Rt a (s a ) 

= 

" llZlZH djpM; s^)Rt b (s b ) + D ab (s> a , s b ) - S^; + m S ds'>D ac {s> a - s>>)Rt b (s b ) 

= t a (s' a )Rt b (s b ) (32) 
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-UlZ'Zi KDM O + dsWUs' a , st)]Rt c (s c ) + V ca (s c ; s' a ) 

= t a {s' a )Rt c {s c ) 

We see that Eq.31 and Eq.32 are, as the sub-section title claims, separable integral 
equations for one-variable, one-parameter functions. But the two alternative forms 
have to be solved for all nine functions and not just for three of them before we 
can even ask the question as to whether they both define the same nine functions 
D xy (s x ; M; s x ). Until we have made this proof, it will be safer to call the solutions 
of Eq. 31 (and their two cyclic permutations) L xy (s x ; M\ s' y ) with the understanding 
that s x is the variable and s' y the parameter while the solutions of Eq. 32 will be 
designated by R xy (s' x ] M; s y ). 

3.3 Solution of the 3-3 problem 

To solve these equations we first define 

r(M-m x ) 2 

t x = ds x t x (s x ; M) (33) 

J(m y +m z ) 2 
AM-m x ) 2 

lxy\Sy) = I ds x L xy (s XJ M , Sy) 

J (m y +m z ) 2 

AM-niyf 

Txyifix) = / dSyR xy (s X j Mj Sy) 

J (m z +m x ) 2 

Then, integrating Eq. 31 over the variables yields 



1 


—t a R 


—t a R 


\ 


( I (s') 


\ 


( 


\ 




t b R 


1 


-i b R 




ha{s'a) 




= t b Rta{s' a ) 




(34) 


t c R 


—t c R 


1 


/ 


V l ca(s' a ) 


J 


V t C Rta(s' a ) 


J 





While integrating Eq. 32 over the variables yields 



1 


—t a R 


—t a R 


( r aa (s' a ) 


\ 


( l 




t b R 


1 


—t b R 


r ab(s' a ) 




= t a {s' a )Rt b 


(35) 


t c R 


—t c R 


1 J 




J 


V ta(s' a )Ri c ) 





These equations are clearly algebraic and, together with the two cyclic permuta- 
tions, define the functions l xy (s y ) and r xy (s' x ) for any value of the entrance or exit 
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parameters allowed by the fixed value of M. Assuming that Eq. 34 has been inverted 
(see below), we then can insert the solution in Eq. 31 to specify the full solution 



L aa (s a ;M; s' a ) 
L ba (s b ; M; s' a ) 
L ca (s c ; M- s' a ) 



M)R[l ba (s' a ) + Us' a )} 
t b (s b ; M)R[t a (s' a ; M) + l aa (s' a ) + l ca {s' a )} 
t c (s c ; M)R[t a (s' a , M) + l aa (s' a ) + l ba (s' a )} 



(36) 



while if Eq. 35 has been inverted we find that 



Raa(s' a , M; s a ) 
R a b(s' a ; M; s b ) 
R ac (s' a ; M; s c ) 



MO + r ac (s' a )]Rt a (s a ; M) 

[t a (s' a ; M) + r aa (0 + r ac (s' a )]Rt b (s b - M) 

[t a (s' a ; M) + r aa (0 + r ab (s' a )]Rt c (s c - M)R 



(37) 



It still remains to prove that, when we complete the system by supplying the remain- 
ing permutations, the two different routes by which we arrive at the solutions define 
the same nine functions. 

One way to do this is to define 



"xy — 



xy 



^xytxRty tyx 

(M-m x ) 2 AM-rriyf 

ds x / ds y L xy (s x ; M; s y ) 

(m y +m z ) 2 J{m z +m x ) 2 
(M-m x ) 2 AM-rriyf 

ds x I dsyR X y(s X j A^, Sy) 

(m y +m z ) 2 J(rn z +rn x ) 2 



(38) 



Then, integrating Eq. 34 over the remaining parameter yields 



-t„R —tnR 



\ 



-t b R 1 — t b R 
-t r R — t r R 1 



(l \ 

vnn. 



ha 

y -l c H -t c H 1 J \ lea / 

while integrating Eq. 35 over the remaining parameter yields 
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We see immediately that the two alternative forms define the same algebraic matrix 



and hence that I 



xy ' yx- 



Because of the symmetry of the driving terms l xy = l yx , 



r xy = r yx and we are free to define 

Zxy = ^xy "^xy lyx ^"yx 

The inversion is algebraic and straightforward fT2|j . Explicitly, with 

R = i 

we find that 



(41) 



(42) 



' z \ 



t a 



1 + t a tb + tbt c + t c t a + 2it a tbt c 



(43) 



Zba itfo tbt c 

v Zca J V ttc t c tb J 

We have now solved the 3-3 problem (with no bound states) because we can 
reconstruct the variable and parameter content of the solution as 

For A = B = C = 

t a (s a , M)6(s a - s' a ) + t a (s a , M)z aa t a (s' a} M) 
t b (s b ,M)z ba t a (s' a ,M) (44) 
t c {s Cl M)z ca t a {s' al M) 

Clearly, the remaining six amplitudes can be written in the same way. 

The Kowalski version of the unitarity proof is now algebraic and trivial, provided 
we use the state normalization which makes two-particle unitarity require that t — t* = 



M aa (s a ; M; s' a ) 
M ba (s b ; M; s' a ) 



M ca (s c ; M; s' a ) 



2itt* and R — R* = 2iRR* , which forces us to choose the constant R 



i. 



3.4 Coalesence, Breakup, Anelastic Scattering coupled to 3 
free-3 free Scattering 

If we now include N a + N b + N c two-body bound states \^ x > with x G a,b,c and 
n x G 1,2, ..,N X the ort honor mality condition must be extended to include 



< s x \/4 > 



OxtJO, 



(45) 
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and the completeness equation becomes 

-{M-m x ) 2 



ds x \s x >< s x \ + >< /CI = 1 (46) 

(m y +m z y 

Although the propagator lZ8 xy remains constant, using our normalizations in this 
extended space, we can add a new dynamical element to the system if we allow (as we 
can, conserving on-shell energy and 3- momentum) the 2-2 scattering operators T a (M), 
whose matrix elements in the 3-3 part of the space are given above, to couple 1,2 states 
(anelastic scattering) to 3 free particle states (breakup) or visa versa (coalesence) via 
the three free particle propagator 1Z preserving probability conservation (unitarity). 
In constructing our specific way of providing this coupling, we have been guided by (a) 
the success of the dispersion-theoretic non-relativistic approximation of single pion 
exchange in predicting the 1 Sq shape parameter in the archetypal strong interaction 
problem (nucleon-nucleon scattering [231 P15|), (b) Faddeev's analysis of the role of the 



"essential singularities" in the non-relativistic case and (c) the success of the "handy 



dandy formula" connecting masses to coupling constants |16|, |22 |. 

We note first that the thresholds for the opening of the elastic and anelastic 
channels can be uniquely specified by the thresholds: 

M nx = m x + aC (47) 

and (assuming no degeneracies and at least one bound state) uniquely ordered by 
some integer parameter 1 < n < N a + iV& + N c . So we can always specify and order 
uniquely the parameters 

M L < M n x , < M Zn <m a + m b + m c - x, x', x" G a, b, c; n < N a + N b + N c (48) 

Consider first only anelastic scattering (M^ < M < m a + mb + m c ). Then, relying 
on Faddeev's insight (b), we will assume that the scattering operator has the matrix 
element 

U£M((*a - 0C) 2 ) < *zr\%\tf >) = (49) 

S xy S nx< iT n x *(k n x *,M)6(M - m x - ^) 
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Here T™ x (k x x , M) is the residue at the pole of the two-body scattering amplitude at 
s x = (/i™ 1 ) 2 and for our scattering length model is given by (cf.Eq.'s 2,3,4,8) 



rs-(A£-, M) = ^[M* - K + ^YW 2 - (m x - (50) 

This is a simple algebraic consequence of the algebraic structure of our two particle 
input model; The ^-function in Eq.49 makes this limit consistent with the orthogo- 
nality (45) and completeness (46) relations. Both below and above breakup threshold 
we model the opening and closing of the bound state vertex by assumung that 

n' 

< ^ x \TZ\ny y > = 5 xy R 
<^ x \K\s y > = S xy R (51) 

n' 

< s x \TZ\fiy y > = 5 xy R 

For M < m a + + m c , taking matrix elements of the Faddeev equations between 
the bound states and inserting the completeness relations as we did in the continuum 
case, we clearly obtain N algebraic equations for N amplitudes where N is the number 
of open channels allowed by the ^-fuctions using the fixed value of M. For one bound 
state in each channel, we obtain a triple of three equations in three unknowns, whose 
solution is of the same algebraic form as those for the z xy exhibited above, with t x 
replaced by r*. These were the equations we had in mind in the first version of 
this paper EH]; the intent of the original paper was to explore breakup and related 



problems in subsequent work(T^]. As already noted, the more complete treatment 
here was needed to meet referee comments, for which we are grateful. 

The postulate of a "channel independent" (indeed constant) "three free particle 
propagator" below breakup threshold is our way of taking over the Faddeev idea that 
"once a scattering occurs, we must allow one of the pair to interact with a third 
particle before anything more can happen" into our on-shell theory. Another way 
to put it is that "once a bound state vertex opens up, another degree of freedom 
must intervene before that pair can interact again." This is familiar in other Faddeev 
contexts; only this particular articulation of the idea is novel. This analysis also 
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shows us how the coalesence and breakup amplitudes couple into the system (cf. 
Eq.'s 49,50,51). 

It is now straightforward to write down the full coupled equations for A, B,C,V 
above breakup threshold for any arbitrary finite number of bound states and reduce 
them to algebraic equations for the constants z X y v ■ Once solved, we can reconstruct 
the full variable content of the amplitudes by writing 



v 

9(M - (m x + fj^))9(M - K + $ )) (52) 
B2l(s x ;M) = t x (s x ,M)zivf 

9(s x - (m a + m b + m c ) 2 ) (53) 
C n x ;{M-s' y ) = T n x *z n x ;t y {s' y ,M) 

9{s' y - (m a + m b + m c ) 2 ) (54) 
V xy (s x ;M;s' y ) = t x (s x ,M)z xy t y (s' y ,M) 

9(s x - (m a + m b + m c ) 2 )9(s' y - (m a + m b + m c ) 2 ) (55) 

Knowing z X y Uv , z X y, z xy and z xy we can now provide predictions for all elastic, re- 
arrangement, three body breakup and 3-3 scattering cross sections for all energy- 
momentum-J = conserving processes over the entire kinematic region M < M t h. 

The two body model used assumes that we know either the binding energy of a 
(single) two body bound state or the scattering length in the isolated two-particle 
system but not both. Our treatment up to this point has implicitly assumed that our 
six parameters m x ,n x , x G a,b,c, are all consistent with known two body data over 
the range of interest. The failure of this connection, either by the two parameters 
(which can be determined by different types of experiment) being inconsistent with 
it, or by the departure of the predicted elastic scattering cross section in the physical 
region from experiment suffices to show that we must enrich the parameter content of 
the model. This would be analagous to what happened historically || in the study of 
neutron-proton scattering and established the spin-dependence of nuclear forces. Pre- 
cise experiments analyzed using related ideas have even demonstrated the existence 



of pions using nucleon-nucleon s-wave experiments below 10 Mev[24, 19|. 
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When it comes to the unique predictions of our three body model, the comparison 
with experiment becomes richer than for the two-body scattering length model. In 
particular we can now predict three scattering lengths or their equivalent and com- 
pare these with low energy s-wave experiments. If any one of these fails to agree with 
experiment, one place to look for an explanation is to postulate a single three-body 
bound state. If our matrix formulation of the Faddeev equations, when analytically 
continued below the threshold for the lowest kinematically allowed 2-2 elastic scatter- 
ing, possesses a homogeneous solution at some value of M, this value of M predicts 
the existence of a 3-body bound state at that invariant energy. As already noted, this 
analytic continuation in our theory is only possible after we embedded our model in 
an appropriate 4-body space. If the state with the right quantum numbers is found, 
but the value of M differs experimentally from that predicted, we can extend our 
model phenomenologically by explicitly introducing an (abc) «-> (abc) channel in ad- 
dition to the three channels we started with. Then we must solve four equations for 
four unknowns, but no new conceptual problems arise. Since direct scattering of three 
free particles to three free particles is usually too difficult to measure, particularly 
at three body breakup threshold, we do not count this as a new piece of available 
experimental information. However, we now have one parameter to explain (if the 
bound state exists at a known mass m a b c ) four experimental numbers, giving a strin- 
gent consistency check even using threshold data. We intend to explore consistency 
conditions on these parameters elsewhere. 

As a significant example, particularly relevant to Tjon's oevre, the n,n and n,p 
singlet scattering lengths and the deuteron binding energy can be used as the three 2- 
2 channel parameters for the n, n,p system. Only two scattering lengths (n, d doublet 
and n, d quartet) are measured in addition to the binding energy of the triton. 4 a n d is 
well predicted while 2 a n d and e t are highly correlated ( "Phillips line" ) for reasons that 
can be understood from a dispersion-theoretic point of view|3|]. Of course Barton and 
Phillips' explanation is consistent with the physics underlying our on-shell approach, 
so we expect our model to achieve a comparable result. In historical fact, it was their 
work which helped us start thinking about the usefulness of a more general approach 
in the first place. 
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It should be obvious from the treatment of the four body problem in our non- 
relativistic paper |H] that the current approach can be readily generalized to relativis- 
tic four particle systems, as we intend to do on another occasion. 



4 SPECULATIONS AND CONCLUSIONS 

The example of the applicability of this model to nuclear physics given in the last 
section hardly begins to suggest the range of problems which we believe can profitably 
be explored using the approach presented here. For instance, the equivalent of the 
relativistic scattering length formula used here was first written down by Bohr in 
1915 16| . Viewed from the current point of view, this makes the hydrogen atom 



a relativistic bound state of a proton and an electron. This suggests looking at the 
three body systems e,e,p (H~), e,p,p (H2) and similar atomic systems using the 
explicit model presented in this paper. 

For strong interactions, we suggest treating the deuteron as a neutron-proton-pion 
bound state. If we include crossing and the Fermi- Yang model for the pion|10j as a 
bound state of a nucleon and an anti-nucleon the usefulness of the approach for deeply 
bound states should become manifest. We suspect that relativistic models of quarks 
and quark confinement could also be attacked using these methods. 

The relativistic scattering length model employed here actually arose in a study 
of the fine structure spectrum of hydrogen ||16|| . In conjunction with combinatorial 
arguments, this model leads to the Sommerfeld formula without any specific use of the 
concept of "spin" , and to an improvement of the lowest order combinatorial calcula- 
tion of the low energy fine structure constant (a e (m 2 e )~ l = 137) by four significant 
figures (to 137.03596...). Similar improvement of our understanding of "bit-string- 
physics" J22] can be anticipated when we make use of the three and four body dy- 



namics adumbered here. We hope that others may be induced to try these simple 
methods and see how far they might lead. 

Because of the occasion to which this issue of Few Body Physics is dedicated, one 



of us (HPN) thinks it appropriate to reiterate here, as was stated long ago [18], that the 



reduction of the three body problem from three to two continuous variables presented 
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by Osborn and Noyes|28[] was first, independently, developed by Ahmadzadeh and 
TjonQ. This is yet another reminder, of which there will be many in this issue, of 
how continuously useful and important Prof. Tjon's dedication to our field has been. 
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